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A study based on ab initio calculations is presented on the structural, elastic, and vibrational 
properties of single-wall carbon nanotubes with different radii and chiralities. These properties 
are obtained using an implementation of pseudopotential-density-functional theory which allows 
calculations on systems with a large number of atoms per cell. Different quantities are monitored 
versus tube radius. The validity of expectations based on graphite is explored down to small radii, 
where some deviations appear related to the curvature-induced rehibridization of the carbon orbitals. 
Young moduli are found to be very similar to graphite and do not exhibit a systematic variation 
with either the radius or the chirality. The Poisson ratio also retains graphitic values except for a 
possible slight reduction for small radii. It shows, however, chirality dependence. The behavior of 
characteristic phonon branches as the breathing mode, twistons, and high-frequency optic modes, is 
also studied, the latter displaying a small chirality dependence at the top of the band. The results 
are compared with the predictions of the simple zone-folding approximation. Except for the known 
defficiencies of the zone- folding procedure in the low- frequency vibrational regions, it offers quite 
accurate results, even for relatively small radii. 

PACS numbers: 63.20Dj; 61.48. +c; 78.30.Na; 71.15.Mb 



I. INTRODUCTION 



Carbon nanotubes have excited a considerable interest 
in the condensed-matter and materials research commu- 
nities in the last few years and much experimental and 
theoretical work has been devoted to them as prototype 
of one-dimensional ordered systems with promising tech- 
nological applicationsa. Electronic transport in conduct- 
ing nanotubes is one of the issues that has attracted more 
attention, especially after the developments that made 
possible the sinthesis of large quantities of_single wall 
nanotubes (SWNT) forming crytalline ropesall. Experi- 
ments showed a peculiar metallic behaviour above 35 K 
which was understood in terms of the coupling between 
the conduction electrons and Jong wavelength twistons, 
i.e., torsional shape vibrationsu. The quantitative under- 
standing and characterization of this and other related 
phenomena require the detailed knowledge of both the 
structure and the vibrations of these tubes. 

The structure of carbon nanotubes is qualitatively well 
known through the simple construction of rolling a per- 
fect graphene sheet, where only one parameter is to be 
determined: the lattice parameter or a bond-length. The 
symmetry of the tubes is less restrictive than in graphite 
and several parameters are needed to determine com- 
pletely the structure. Among other things, these param- 
eters define the differences between inequivalent bonds, 
which condition the position of the Fermi surface in the 
conducting armchair tubes. It is very difficult to ob- 



tain direct experimental information for the structure, 
and j*ry little theoretical information has been given so 
faritl. 

Besides possible nanotechnological applications, car- 
bon nanotubes are promising candidates for composite 
materials where their low weight and very high Young 
modulus can be of use. Their clastic properties have thus 
received considerable attention as well. The nanotube's 
unusual strength arises from a combination of high stiff- 
nesS|-and a extraordinary flexibility and resistance to frac- 
tureErtL2l. The Young modulus of SWNT has been mea^. 
sured from the amplitude of their thermal vibrationsliiJ 
and by measuring the bending fpsce of a pinned nanotubc 
by an atomic force microscope£3. The experimental re- 
sults give a Young modulus in the range of one TPa, 
similar to the one of graphite when pulled parallel to the 
sheets, but the experimental uncertainity is WMjtfii high. 
There is also a dispersion of theoretical valuesfflll3il3 in 
the literature for this quantity corresponding to different 
approximations but also to different definitions of the ef- 
fective sectional area. 

Raman experiments^ offer valuable information for 
the vibrations of ropes of single-wall armchaii (n, n) 
tubes. The use of resonant Raman scattering^ allows 
to discriminate the vibrations stamming from tubes of 
different diameters by looking mainly at the Ai ff breath- 
ing mode that exhibit a strong dependence on tube- 
diameter. Furthermore, the optical Fi2 g phonon peak 
shows size-dependent multiple splittings that are nearly 
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independent of chkality and can be used to determine 
the tube-diameterE3. Raman scattering, however, is lim- 
ited to the neighborhood of the center of the Brillouin 
zone, and additional selection rules limit the number of 
vibrations detectable to seven per tube type. The theo- 
retical characterization of the vibrational modes can com- 
plement this knowledge allowing to correlate tube diam- 
eter to a specific vibrational frequency. There is theoret- 
ical information available based on empirical force con- 
stantslHl and tight-binding HamiltoniansEjO, but a first- 
principles theoretical reference is still lacking£3 (except 
for the breathing modeQ). 

The strong similarity of the chemistry of carbon nan- 
otubes to graphite allows theoretical analyses to be done 
based on empirical methodologies imported from studies, 
on graphite. They range from the direct zone-foldingH 
of the results for graphite to the quantum-mechanical 
studies based on tiabipbinding Hamiltonians fitted to 
graphite propertiesE3Eilc3. Effective interatomic poten- 
tials^, iarce-constant modelsEj or non-orthogonal tight- 
bindingt£l have also been tried. The performance of the 
different techniques varies, from the qualitative picture 
offered by zone folding, with intrinsic defficiencies for 
low frequencies, to the very quantitative results of tight- 
binding approaches. The curvature of the tubes, how- 
ever, disturbs the chemistry in a way that can cause 
the deviation, from the graphite based description, for 
narrow tubes. Zone folding and force-constants neglect 
curvature all together. Model potentials can only ac- 
count for the different distances among the atoms. Tight- 
binding captures part of the chemical strain through the 
geometry-dependence of its electronic matrix elements, 
even though their absolute value depend on the electronic 
structure of graphite. It is then important to be able to 
compute the different properties for any tube radius using 
a tool that does not depend on a fit to graphitic proper- 
ties, so as to study the narrow-tube properties with the 
same degree of accuracy as the wide ones. 

In the present work we present an ab initio study of the 
structural, elastic, and vibrational properties of single- 
wall carbon nanotubes for different diameters and chiral- 
ities to address the points mentioned above. The behav- 
ior of the different properties is monitored versus tube 
radius, and the performance of the different approxima- 
tions used in the literature is studied to ascertain on their 
different ranges of applicability. Section II describes the 
methodological details. The results and discussion are 
presented in Section III, to finish with the concluding 
remarks in Section IV. 



II. COMPUTATIONAL SCHEME 

The first-principles scheme used in this study was de- 
scribed in detail elsewhereE3, where it was already ap- 
plied and tested on large fullerene molecules. It is based 
on the Local Density Approximation (LDA) to Density 



Functional TheoryEa. Core electrons are replaced by non- 
local, norm-conserving pseudopotentialsEH, whereas va- 
lence electrons are described within the Linear Combi- 
nation of Atomic Orbitals (LCAO) approximation. In 
this work we have used a minimal basis set of one s and 
three p orbitals per carbon atom. Although this basis is 
certainly not complete, it provides a sufficiently accurate 
description of the systems and effects that we intend to 
study. 

The radial shape of the atomic basis functions is chosen 
according to the prescription of Sankey and NiklewskiEa. 
This consists in using the solution of the atom with the 
pseudopotential as a basis for the LCAO calculation; 
these pseudoatomic orbitals (PAO) are calculated with 
the boundary condition that they vanish outside a given 
radius r c . The PAO's are therefore slightly excited, since 
their kinetic energy is increased due to the vanishing 
boundary condition. They are strictly localized, which 
is very advantageous from the computational point of 
view. For a review of the quality of these PAO bases for 
electronic structure calculations, we refer the reader to 
Ref. In this work, we have used a radial cutoff of 
r c = 4.1 Bohr. 

The use of atomic orbitals, and a combination of effi- 
cient techniques, allow us to calculate the LDA Hamil- 
tonian with an order- N effort (i.e., a cost that scales 
linearly with the number of atoms, both in time and in 
memory)E3. This makes possible to reach system sizes 
with a much larger number of atoms than the standard 
techniques. Some of the Hamiltonian matrix elements 
are computed by interpolation of pre-calculated two and 
three center integral tables, whereas others are obtained 
by direct integration in a real space grid. The fineness of 
the grid is expressed by the maximum kinetic energy of 
a plane wave that can be represented in it, as is usually 
done in plane wave calculations. Note that, in our case, 
the grid is used to represent the charge density, and not 
the wavefunctions. In this work we have used a cutoff of 
60 Ry for the grid integrations. 

The solution of the Hamiltonian matrix can be done by 
straight diagonalization (which is an 0(N 3 ) operatkjpA, 
or using recently developed order-iV techniquesEIEil, 
which are very advantageous for systems with large num- 
bers of atoms (tipically, the crossover between the two 
solutions is around 100 atoms, for calculations with min- 
imal bases) . The systems under study in this work are in 
the range between 80 and 200 atoms, and therefore diago- 
nalization of the Hamiltonian is feasible and competitive 
with an order-iV technique. We have used diagonaliza- 
tion throughout this work. 

Calculations were performed for the following tubes: 
(4,4), (6,6), (8,8), (10,10), (10,0), and (8,4). They were 
all considered as isolated, infinitely long tubes. For that 
purpose, we used periodic boundary conditions on a su- 
percell geometry with sufficient lateral separation among 
neighboring tubes. For the purpose of sampling the Bril- 
louin zone in the direction of the tube axis, as well as for 
the computation of the force constant matrix for phonon 



2 



calculations (see below), we used supercells containing 
several unit cells in the axis direction. The supercells 
consisted of five unit cells for the armchair (n, n) tubes 
(80, 120, 160, and 200 atoms for the (4,4), (6,6), (8,8), 
and (10,10) tubes, respectively), three unit cells for the 
zig-zag (10,0) tube (120 atoms), and one unit cell for the 
chiral (8,4) tube (112 atoms). The supercell length in the 
tube-axis direction is similar for all of them, the number 
of atoms changing because of the different tube radii. 
Only the T point of the supercell was used in the phonon 
calculations, although more complete fc-point samplings 
were tested showing no significant differences with the F 
point results. 

The atomic structures of the tubes were obtained by 
careful minimization of the total energy by means of the 
calculated Hellman-Feynman forces, including Pulay-like 
corrections to account for the fact that the basis set is 
not complete and moves with the atomso. The mini- 
mizations were performed using a dissipative molecular- 
dynamics algorithm which allows the geometry optimiza- 
tion with no symmetry conditions imposed. The resid- 
ual forces in the optimizations were always smaller than 
0.04 eV/A. The energy was minimized and the structure 
relaxed for different values of the lattice constant along 
the tube axis. This procedure allows us to determine the 
most stable lattice constant along the tube, as well as the 
the Young modulus and the Poisson ratio. 

For the phonon calculations, we first obtain the force 
constant -matrix in real space using a finite difference 
approachBEi We used finite atomic displacements of 
0.02 A, and the force constants were taken as the aver- 
age of the results obtained with positive and negative dis- 
placements, to elliminate anharmonic effects. All atoms 
are equivalent by symmetry, in the nanotubes, and there- 
fore we only calculated the force constants for one of the 
atoms in the supercell, and generated the rest of the ma- 
trix using the symmetry operations. 

The force-constant matrix has to be computed between 
a given atom and all the rest in the system. However, it is 
known that the force-constants decrease rapidly with dis- 
tance (in non-polar systems), so that only the elements of 
atoms sufficiently close need to be computed. To do so, 
we set up a supercell large enough that a sufficient num- 
ber of neighbors in the tube axis direction is included. 
It must be also kept in mind that, in the supercell ge- 
ometry, a given atomic displacement in the central cell is 
always accompained by the same displacement of all the 
images. The supercell must therefore be large enough so 
that the effect of the image displacement is neglegible. 
The cells discussed above were built with these obser- 
vations in mind. Once the force constant matrix in real 
space has been obtained, we calculate the dynamical ma- 
trix in reciprocal space, and diagonalize it to obtain the 
phonon modes and frequencies as a function of the 1-D 
crystalline momentum vector in the direction of the tube. 



III. RESULTS AND DISCUSSION 
A. Structural properties 

We first study the equilibrium structural properties 
of the nanotubes. As a reference, we have first com- 
puted the equilibrium bond distance for a single graphene 
plane, for which we obtain a value of d= 1.436 A, close to 
the experiment alti-l (1.419 A) and LDA plane- wave cal- 
culationsE3 (1.415 A) in graphite. Our value is slightly 
larger than the one obtained with plane- waves, the small 
difference being due to the basis set. For all the tubes 
studied here we have found that the average carbon bond- 
length is within 1% of the graphitic value. No apprecia- 
ble symmetry-breaking distortions have been observed. 
However, by symmetry, for the (n, n) and (n, 0) tubes 
there are two inequivalent bonds, three in general, for 
chiral tubes. Fig. [l] (a) shows the value of these two 
bond lengths as a function of tube radius (in units of the 
graphene bond length) for the (n, n) tubes studied here. 
The differences between these bonds and graphite, and 
between the two kinds of bonds, are small but signifi- 
cant. Two effects are apparent: (i) both bond-lengths 
increase (as compared with the graphene reference) with 
decreasing tube radius; (ii) the difference between them 
also increases with decreasing radius. Both effects can 
be easily understood in terms of the rehibridization and 
the weakening of the it bonds induced by the curvaturdj. 
As shown in Fig. |l| (a), the longer bond is the one per- 
pendicular to the tube axis in the (n,n) tubes. This is in 
contrast with the results of Ref. |7|. This discrepancy, is 
due to the small k-point sampling used in that workEx 

Fig. [l] (b) shows the variation of the bond angles with 
the tube radius for the (n, n) tubes. As for the bond 
distances, there are two inequivalent bond angles in (n, n) 
tubes. The behavior of the bond angles is very similar 
to the one expected from ideal 'rolling' of a graphene 
plane to form the tube, also shown in Fig. [I] (b). In this 
case, one of the angles would maintain a value of 120°, 
whereas the other would decrease for smaller tube radii, 
leading to the increasing tube curvature. In the fully 
relaxed structures, both angles are smaller than 120°, so 
that the curvature stress is more distributed around the 
tube area. Since the ideal 120° angle now shares some 
of the cost of the curvature, the other angle increases 
slightly from its ideal value. Nevertheless, we see that 
deviation from the ideal behavior is only relevant for the 
tubes with small radius, being almost negligible for tubes 
larger than (6, 6). 

As a consequence of the increase in bonding distances 
shown in Fig. [l] (a), the actual tube radii are slightly 
larger than those resulting from an ideal rolling of the 
graphene sheet. This is plotted in Fig. |l| (c). Again, this 
deviation is more pronounced for the tubes with smaller 
radius, and the behavior approaches the ideal for large 
radii tubes. 

The structural distorsions affect the electronic struc- 
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FIG. 1. Structure variations with tube radius, (a) Lengths 
of the two inequivalent bonds in (n, n) tubes in units of the 
bond length in graphite. Circles stand for the bond perpen- 
dicular to the tube axis, and the squares for the other inequiv- 
alent bond-length, (b) The two inequivalent bond angles in 
(n, n) tubes. Continuous lines show the angles resulting from 
an ideal graphene rolling, (c) Tube radius dilation as com- 
pared to the ideal radius, (d) Position of the Fermi surface 
(kp) in the Brillouin zone. Squares and circles stand for the 
ideal and relaxed structures, respectively. The continuous line 
indicates the ideal zone-folded value. 



ture of these metallic tubes. Although the study of the 
electronic properties is not the focus of this work, it is 
worth discussing what is the effect of the structural pa- 
rameters on the Fermi level, since this has implications 
for the physics of these systems, for example, in the inter- 
pretation of very recent scanning tunneling spectroscopy 
experiments on short SWNTEZI. In a zone-folding model 
based on the graphene band structure, the Fermi level 
of all (n, n) tubes would be located at a wave vector of 
|f (with L being the lattice constant in the tube axis). 
For the real tube structure, the position of the Fermi 
level will change due to the reduced symmetry, with its 
two inequivalent bonds (those perpendicular to the tube 
axis, and those non-perpendicular). In a simple Hiickel 
model this reduced symmetry is described with two dif- 
ferent hopping interactions t± and t' , both different from 
the ideal graphene hopping t. The Fermi wave vector in 
this case is &p = %■ arccos(|p-). In real tubes, the Fermi 
surface moves because of two reasons: (i) The rolling of 
the graphene plane to form the tube originates a change 
of the electronic potential in the inequivalent bonds of 
the tube, even in the case where the structure is taken 
as the ideal graphene rolling, i.e. with both bonds kept 
at the same graphene bond length, {it) The difference 
in bond lengths for the two inequivalent bonds in the re- 
laxed structure change the hopping matrix elements. In 
order to distinguish these two effects, we show in Fig. |l| 
(d) the calculated values of the Fermi wave vector for 
the fully relaxed structures, and for the ideal graphene 
rolling tubes, obtained with our LDA-LCAO formulation. 
We see that the main effect of the Fermi level shift is 
the rolling to form the tube, whereas the difference in 
bond lengths brought by the structure relaxation gives 
only a minor correction tending to bring kp closer to 
the graphene value. We see that the deviation from the 
graphene Fermi level are relatively large, increasing with 
decreasing tube radii. 



B. Elastic properties 

1. Strain energy 

Fig. |^ (a) shows the strain energy per atom (energy rel- 
ative to a planar graphene sheet) as a function of the ra- 
dius of the tube. The data follows quite well the behavior 
expected from classical elasticity theoryEj, E st — C/r 2 , 
where r is the radius of the tube and C is a constant 
that depends on the Young modulus Y and thickness h 
of the wall in a model tube: C — Yh 3 a/24. A least 
squares fit to the results of the (n, n) tubes yields a value 
of C = 2.00 eVA 2 /atora., in very good agreement with 
recent LDA calculations!! For the other two tubes stud- 
ied, (8,4) and (10,0), we obtain slightly larger values (2.15 
and 2.16 eVA 2 /atom, respectively). 

It is to some extent suprising that the predictions from 
elasticity theory are so closely followed by the detailed 
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ab initio calculations. In fact, the assignment of a thick- 
ness h for a single atomic layer is not a well defined proce- 
dure. Adams et alSB provided an alternative explanation 
based on microscopic arguments. They use a very sim- 
plified model in which the energetics of many different 
fullerene structures depend on a single structural param- 
eter: the planarity <f>^, which is the angle formed by the 
7r orbitals of neighbour atoms. Asuming that the change 
in total energy is mainly due to the change in the t v in- 
teraction between these orbitals, and that this change is 
proportional to cos0„., the r~ 2 behaviour is predicted. 
Adams et al. obtained a value of C — 2.12 eVA 2 /atom 
using non-selfconsistent first- principles calculations. Pre- 
vious calculations using Tersoff and Tersoff-Brenner po- 
tentialsE3 predict the same dependence and gave a value 
of C w 1.5 eVA 2 /atom and C ~ 1.2 eVA 2 /atom, respec- 
tively. 

We note in Fig. || (a) that the (n, n) tubes are energet- 
ically more stable as compared to other chiralities with 
the same radius. This difference is, however, very small 
and will decrease as the tube diameter increases due to 
the r -2 behaviour. This is expected, since in the limit of 
large radii the same graphene limit is obtained, regardless 
of chirality. 



2. Young modulus 

Measurements of tha-amplitude of the intrinsic thermal 
vibrations of the tubestil allow an undirect determination 
of the Young modulus, yielding an average value of 1.8 
TPa. There is a large uncertainty in the results, however, 
with values ranging from 0.40 TPa to 4.15 TPa. Further- 
more, measurements of the restoring force on bent nan- 
otubes with an atomic force microscope give an average 
value of 1.28±0.59 TPa. Since the corresponding elastic 
constant for graphite (en) is also of the order of 1 TPa, 
it is still a question whether the rolling of the graphene 
plane to obtain nanotubes increases or decreases its stiff- 
ness. 

Theoretical calculations can help in solving the issue, 
but the calculations performed so far also show a great 
dispersion in their predictions. Some theoretical esti- 
mates report values of the order of 5 TPao based on an 
empirical Keating force-constant model for finite-capped 
(5, 5) tube. This unreasonably high value can be due 
to the small size of the aggregates used to describe, the 
tube, with only up to 400 atoms. Yakobson et alB, in 
a study of the structural instabilities of SWNT for large 
deformations, and using Tersoff-Brenner potentials, ob- 
tain an estimate for the Young modulus of about 5.5 
TPa, by fitting their results to the continvmpi elasticity 
theory. However, D. H. Robertson et alE3, using the 
same potential, report a small weakening of the stiffness 
of nanotubes as the diameter decreases, and therefore a 
Young modulus for the nanotubes smaller than the one 
of graphene. They find a systematic dependence with 
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FIG. 2. (a) Strain energy versus tube radius. Solid line 
drawn across the (n, n) data corresponds to a least square 
fit to the C/r 2 behaviour. The two C/r 2 functions passing 
through the (8,4) and (10,0) data are also shown (in this scale 
they appear as one line). The r~ a behaviour is clearly shown 
in the inset. The value obtained for a from the logarithmic fit 
is 2.05 ± 0.02. (b) Young modulus versus tube radius. Open 
symbols for the multi-wall geometry, and solid symbols for 
the single-wall-nanotube crystalline-rope configuration. The 
experimental value of the en elastic constant of graphite is 
also shown, (c) Poisson ratio versus tube radius. 
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the chirality which, although beintsmall, increases with 
decreasing tube radius. J. P. LuEJ, using a force con- 
stant model fitted to reproduce the phonons and elastic 
constants of graphite, obtains elastic properties which 
are essentially independent of helicity and tube radius, 
and comparable to those of graphene (with values of 
the Young modulus below 1 TPa). Finally, recent tight - 
bindingt£l calculations give values of the order of 1 TPa 
for the SWNT, quite insensitive to the chirality of the 
nanotube, and mainly determined by the tubule diame- 
ter, approaching the graphitic limit for diameters of ~ 
1.2 nm. 

Part of the discrepancies in the theoretical results dis- 
cused above is merely due to a different definition of the 
Young modulus in these systems. From the point of view 
of elasticity theory, the definition of the Young modulus 
involves the specification of the value of the thickness of 
the tube wall. As discussed in the previous Section, it is 
not clear how to define this width for a SWNT, where the 
wall is composed of only one shell of atoms. The anoma- 
lously large value obtained by Yakobson et aln is due to 
an assignment of a value of h = 0.6 A for the thickness of 
the grapheme plane, which is obviously too small. Other 
authorsr 5 H 16 l have used the graphite interlayer spacing of 
3.4 A in their calculations. 

In order to avoid this definition problem, we have anal- 
ized our results on the elastic stiffness of the nanotubes 
using the second derivative of the strain energy with re- 
spect to the axial strain: To obtain this quantity 
for the different tubes, we have performed structural re- 
laxations for the nanotubes, subject to deformations be- 
tween -1.0% and 1.0%, at intervals of 0.13%. The total 
energies for deformations in the interval (-0.75%, 0.75%) 
where fitted to a third order polynomial, and was 
obtained from the second derivative at zero strain. The 
results are shown in Table ||. We have estimated the nu- 
merical error of these results by taking different intervals 
in the fitting procedure: (-0.5%,0.5%) and (-1.0%, 1.0%), 
and obtain that the typical uncertainty of the calculation 
of ^ is of the order of 10%. From Table Q we see that 
the average value in the tubes is about 56 eV. The vari- 
ation between tubes with different radii and chirality is 
very small, and always within the limit of accuracy of the 
calculation. We therefore can conclude that the effect of 
curvature and chirality on the elastic properties is small. 

For comparison, we have also calculated the values of 
for a single graphene sheet. Three calculations were 
done, with different supercells. In all the cases, only the 
r point of the supercells was used reproducing the k- 
point sampling of the tube calculations. Two of them 
were the unrolled equivalents of the supercells used for 
the (4,4) and (10,10) tubes. These cells were used to al- 
low for a direct check of the curvature effects, comparing 
the planar and tubular geometry. The third one was a 
rectangular (nearly square) cell of 17.41 Ax 17.23 A, with 
112 atoms, thus providing a more uniform fc-point sam- 
pling, and, therefore, a more confident estimation of the 



TABLE I. Calculated values of d 2 E/d 2 e for different tubes. 
The values for a graphene plane obtained with different su- 
percells are also shown (see text for details). 
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elastic constant. The results (also displayed in Table |) 
clearly show that there are no appreciable differences be- 
tween the results obtained for the nanotubes and those 
of graphene, the differences being within the uncertainty 
of the calculation. These data confirm that the effect of 
curvature on the Young modulus of the SWNT is small 
down to radii of the order of, at least, 2.8 A. 

Our results are in .good agreement with those obtained 
by Robertson et al.E3 using Tersoff-Brenner potentials, 
who find values around 59 eV/atom, with very little 
dependence with radius and/or chirality. Futhermore, 
we can obtain an experimental estimate of this quan- 
tity using the elastic constant^ en = 1.06 TPa of bulk 
graphite, from which we obtain ~ cnV a = 58.2 
eV/atom (where V a is the atomic volumen in graphite). 
This value agrees well with th&*esults obtained here, and 
with those of Robertson et a/B3 

If one insist in calculating the Young modulus using the 
standard bulk definition, instead of the well defined sec- 
ond derivative used above, one must choose a definition 
of the effective area per carbon atom. Here we calculate 
the Young moduli considering two different geometries: 
(i) a multiwall like geometry, in which the normal area 
is calculated using the wall-wall distance as the one in 
multiwall tubes, which is very approximately equal to 
the one of graphite, i.e. the effective thickness of the 
the tube wall is taken to be 3.4 A, and (ii) a rope con- 
figuration of single wall tubes, where the tubes would 
be arranged forming an hexagonal closed packed lattice, 
with a lattice constant of (2r + 3.4 A), being r the tube 
radius. The results are shown in Fig. g (b). For the crys- 
talline rope geometry, the decrease of the Young modulus 
with increasing the tube radius is due to the quadratric 
increase of the effective area in this configuration, while 
the number of atoms increases only linearly with the tube 
diameter. The computed values for the SWNT ropes are, 
however, still very high, even comparing with other car- 
bon fibers.cil 
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FIG. 3. (a) Calculated phonon structure for a graphene sheet, (b) Zone folding result for the (4,4) tube, obtained from the 
graphene ab initio phonons in panel (a). Ab initio dispersion relations for the (4,4) (c), and the (10,10) (d) nanotubes. In (c) 
and (d), thicker lines are used to mark two special branches: the acoustic band is a twiston mode (torsional shape vibrations), 
the other (with finite frequency at F) is the breathing mode. 



3. Poisson ratio 



C. Vibrations 



The Poisson ratio v is given by the variation of the 
radius of the SWNT resulting from longitudinal defor- 
mations along the tube axis: 



Ar 

r 



Al 

'T 



(i) 



where I is the tube length. We have calculated v for the 
tubes under study, and find that in all cases the Poisson 
ratio is positive: an elongation of the tube reduces its 
diameter. The results are shown in Fig || (c). We obtain 
values around v = 0.14 (from 0.12 to 0.16) for the arm- 
chair (n,n) tubes, and a little larger for other chiralities: 
0.19 for (10,0) and 0.18 for (8,4). The uncertainty of the 
obtained values is of the order of 10%. These results re- 
veal a slight decrease of the Poisson ratio with the tube 
radius, and a stronger dependence with chirality. Our 
results are close, to the value of v = 0.19 obtained by 
Yakobson et alB using Tersoff-Brenner potentials, but 
considerably smaller than the value v = 0.28 given by 
LuE3 with a force constant-model and v = 0.26 from 
a tight-binding calculation. lj The corresponding r pagn i- 
tude along the basal plane in graphite is v = 0.16.e3Ej 



In our study of the vibrational properties of SWNT, 
we have first computed the phonon spectrum of a single 
graphene plane, which is shown in Fig. |^ (a). This will 
serve as a test of the accuracy of the calculation method, 
and as a reference for the interpretation of the nanotube 
results. Also, the graphene phonon structure is needed 
to obtain nanotube phonons within the zone-folding ap- 
proach. 

The calculation of the graphene phonons has been per- 
formed using the nearly square supercell of 112 atoms 
described in previous sections. The calculated ab initio 
phonon dispersion curves are, in generaL-jn quite good 
agreement with experiments for graphitecj The most re- 
markable disagreement is for the frequency of the higher 
optical bands, which at T is around 1690 cm" 1 in our 
calculations, and 1580 cm" 1 in the experiment. This dif- 
ference is attributed to the use of a minimal sp 3 basis 
in the calculation. However, our calculations reproduce 
well the overbending of the LO band (i.e., the highest fre- 
quency is not at the T point, but at intermediate points 
between FM and TK). This has important consequences 
in the phonon structure of the nanotubes, especially in 
the zone-folding scheme, where the overbending leads to 
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Raman active modes in the nanotubes which are higher 
in frequency than those of graphene. Also, our results 
show good quantitative agreement with experiment for 
the lower frequency band. For example, for the out of 
plane transversal optical mode (labeled in the Fig. || (a) 
as ZO), we obtain a frequency of 861 cm -1 , to be com- 
pared with the corresponding infrared actived mode in 
bulk graphite of 868 cm" 1 . Acustic bands are also well 
reproduced. 

The main features in the phonon dispersion of Fig. || 
(a) are in gpoeLagreement with other density functional 
calculations 4 ^r^, but differ from the phonon .bands ob- 
tained with empirical force constants models ,EjeZI espe- 
cially around the M point. 

Sound velocities are extracted from the slope of the 
acoustic branches. We obtain 24 km/s and 18 km/s for 
the LA and the in plane TA branches, respectively. These 
results agree very well with the sound velocities that can 
be extracted from the experimental data of Ref. [i~J for 
the graphite (0001) surface phonons, rj 24 km/s and 
« 14 km/s. The out of plane transversal band (ZA in 
the Fig. U (a)) has a zero sound velocity. Fitting the fre- 
quencies below 80 cm -1 to a parabolic function oj = Sq 2 , 
a value of i5 » 6 x 10 -7 m 2 s~ 1 is obtained, in very good 
agreement with previous estimations^. The sound veloc- 
ities of the LA and TA branches allow us to calculate back 
the in-plane stiffness (71.6 eV/atom) and the shear mod- 
ulus (40.3 eV/atom) of the graphene sheet, respectively. 
Both values are higher than the ones obtained directly 
(see previous subsection) due to error propagation: the 
elastic constants depend quadratically on the sound ve- 
locities, and these are delicate to obtain due to numerical 
problems close to the T point. The other way around, the 
sound velocity calculated using the experimental basal 
plane shear modulus of bulk graphite (c66 = 0.44 TPao) 
is 14 km/s, which is indeed quite close to our value. 

From the quadratic behaviour of the ZA band we can 
also estimate the energy necessary to roll up the graphene 
plane to form the tubes. It can be easily shown that the 
strain energy per carbon atom can be approximated by 
E s t = C/r 2 where C = S 2 m c /2, being m c the carbon 
atom mass and r the tube radius. The value obtained 
for C is 2.3 eVA 2 /atom, which is only a 7% higher than 
the calculated constant for the (10,0) and (8,4) tubes, 
and 13% higher than the same constant for (n, n) tubes. 

We next investigate the phonon structures of the nan- 
otubes considered here. Let us first comment on the nu- 
merical errors in the calculated vibrational frequencies. 
These are mainly originated from the residual forces in 
the structural relaxation, the finite atomic displacements 
in the force constant matrix calculation, and the finite 
grid utilized in the integration of the Hamiltonian ma- 
trix elements. We have estimated this numerical error 
from the differences in the frequencies obtained for the 
(10,0) tube by displacing two different (but equivalent by 
symmetry) atoms for the calculation of the force constant 
matrix. We find that the error is about 10 cm -1 for the 
high part of the spectrum (frequencies higher than 1300 
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FIG. 4. Ab initio frequencies of the Raman and infrared 
active modes for a series of (n, n) nanotubes, and for the 
(10, 0) nanotube. Filled symbols indicate the breathing mode. 



cm -1 ), and could increase up to 30 cm -1 for some of 
the lower branches. The breathing mode and the acus- 
tic bands are more stable, showing errors of about 15 
cm . Uncertainties of the same magnitude have been 
reported by other autliops, performing similar calculations 
on graphitic systemsJlj'H 

Fig. H (c) and (d) show the calculated ab initio ID 
dispersion relations for the armchair (4,4) and (10,10) 
tubes. The zone-folding results for the tube (4,4), ob- 
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tained from the graphene ab initio dispersion relations, 
are also displayed in the panel (b) of the same figure. 
The difference between the ab initio and the zone-folding 
frequencies are a consequence of curvature and relaxation 
effects, and are therefore a measure of the their impor- 
tance in the phonon spectrum. 

From the results of Fig. |] we see that, apart from some 
small differences which we will analyze in the following, 
the general agreement between the ab initio results and 
the zone-folding predictions is considerably good. This is 
the case even for the tubes with smaller radii, where the 
curvature effects in the phonon frequencies are expected 
to be more important, and the zone- folding scheme could 
start to break down. The agreement is particularly good 
for the upper part of the spectrum, and worsens for de- 
creasing frequencies. This is mainly due to the failure 
of the zone-folding approach to describe the breathing 
modes and two of the acoustic bands of the tube (those 
corresponding to motion in the directions perpendicular 
to the tube axis). In particular, within the zone- folding 
scheme, the breathing mode appears with zero frequency, 
and the two translational modes appear with finite fre- 
quency. It should be noticed that these deficiencies can 
be correctedEij, and an analytical expresion can be ob- 
tained for the breathing mode frequencies making similar 
assumptions as those made in the zone-folding scheme: 
the_use of force constants from a graphene plane. Jishi et 
alx3 showed that this model predicts a 1/r dependence 
of the frequency of this mode, regardless of chirality. 

The deficiencies mentioned above are absent in the 
ab initio results, without having to resort to additional 
corrections. In Fig. || (c) and (d) we display with thicker 
lines the phonon bands associated with the breathing and 
twiston modes for the ab initio results. 

The twistons are torsional acoustic modes which have 
been proposed to be of relevance for the peculiar lin- 
ear dependence of the of the electrical resistivity with 
temperature in the metallic (n, n) tubes.u These vibra- 
tional modes break the reflexion symmetry of the tubes, 
and open a gap at the Fermi level, producing a strong 
electron-phonon coupling, key to the understanding of 
that behavior. The sound velocity of these modes is 
therefore important for the electronic properties of the 
(n, n) tubes. Our results indicate that the twiston mode 
sound velocity is lower than the corresponding value ob- 
tained for graphene (TA band) for all the studied tubes, 
and slowly diminish with decreasing tube radius. For the 
(10,10) tube, the twiston sound velocity is 15 km/s, i.e. 
a 17% lower than the value found in the graphene plane, 
and for the narrower tube (4,4) the value decreases to 
13 km/s. 

Fig. |J shows the Raman and infrared active modes at 
the T point .j-These have been assigned according to the 
D n d groups,E3 which predict 7 infrared active modes and 
15 Raman active modes for the (n, n) tubes, and 8 in- 
frared and 16 Raman actives modes for the [n, 0) regard- 
less of their radius. We note that, for the (10,0) tube, the 
highest Raman active mode in the ab initio calculation 
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FIG. 5. Ab initio breathing mode frequencies as a function 
of the inverse tube radius, for (8,4), (10,0) and five (n, n) 
tubes. The continuous line is a linear fit to the data ex- 
cluding the (4,4) and (5,5) tubes. Bpt-dashed line shows the 
behaviour obtained by Jishi et alEB using a force-constants 
model. Dashed and dotted Ames show the result of the LDA 
calculations by Kiirti et al.u for the (n, n) and (n, 0) tubes, 
respectively. In this scale the dotted line is hidden by the 
continuous one. Error bars as estimated (see the text). 



has a frequency larger than the correspoding frequencies 
for the (n, n) tubes. In our calculation this highest Ra- 
man active mode has a frequency even larger than the 
predicted by the zone-folding scheme. This is in contrast 
with the findings for the (n, n) tubes. In the (10,0) tube, 
this mode corresponts to an optical vibration in which the 
displacement vector is parallel to the tube axis. Although 
the position of this mode is quite sensitive to the numer- 
ical precision of the calculation (and in particular to the 
Brillouin zone sampling utilized), we can conclude that 
this tube presents optical frequencies which are higher 
than those of (n, n) tubes with comparable radii. 

We have also analyzed the dependence of the breath- 
ing mode frequency with the tube radius and chirality. 
The results are shown in Fig. ^[ The breathing mode is 
a A\ g symmetry mode (being therefore Raman active), 
in which there is a monopolar inward and outward vi- 
bration of the atoms. A simple approach, based on the 
force constants derived from the graphene plane, predicts 
a change of the frequency of this mode as A/r, inde- 
pendently of chirality, where r is the tube radius and 
A=1092 cm-T 1 A (Ref. ||). Recent LDA frozen phonon 
calculations!! of this mode confirm the prediction of the 
r _1 behaviour, the constant A having a weak dependence 
on chirality, 1180 cm" 1 A and 1160 cm -1 A for (n, n) and 
(n, 0) , respectively. Our calculations confirm the previ- 
ous results, indicating once again that the effect of cur- 
vature on the value of the force constants is small, even 
for the small radii tubes considered here. Only the (4,4) 
tube presents an important deviation from the predicted 
behaviour, with an appreciable decrease in the breath- 
ing mode frequency. This effect is already noticeable in 



9 



■r- 

O 
P 




Par. axis 



( c ) Radial Tangential Par. axis 




500 1000 1500 

Frequency (cm -1 ) 



FIG. 6. Vibrational density of states for (a) graphene, (b) 
the (4,4) tube in the zone folding approximation, and (c) the 
ab initio results for the same tube. The curves are decom- 
posed in the different directions of the displacement vector. 
For graphene, Z indicates modes perpendicular to the plane, 
and X and Y within the plane. For the (4,4) tube, the modes 
are decomposed in radial, tangential and parallel to the tube 
axis. 



the (5,5) tube, although to a smaller extent. This effect 
can be understood as a consequence of the hibridization 
changes and the decrease of the n interaction, induced 
by the curvature. Ar 1 fit of the results for tubes with 
radius greater than 3.8 A gives a value of 1160 cm -1 A 
for the constant A, in very good agreement with Ref. [?], 
as clearly shown in Fig. || The possible chirality depen- 
dence of the breathing mode, if any, is well below the 
resolution of our data. As pointed out in Ref. [t], the 
value of A can be estimated from the stretching constant 
of the graphene plane neglecting all the possible effects 
of curvature. Taking our calculated value of 60 eV/atom 
(see previous sections) for this elastic constant, we obtain 
a value A=1166 cm -1 A. 

One of the most important differences between the 
ab initio and zone-folding frequencies is a general soften- 
ing of frequencies when curvature effects are taken into 
account. This is specially clear for the higher frequency 
bands, but is also observable for the intermediate fre- 
quencies. The shift is not uniform, and therefore the 
distribution and, in some cases, the ordering of the Ra- 
man and Ir active modes are affected by curvature effects. 
This effect is most evident in the (4,4) tube, where the 
curvature is more pronounced. We show in Fig. ^| the 
vibrational density of states for the (4,4) tube, compar- 
ing the ab initio and the zone-folding results. The figure 
shows the decomposition of the modes in radial and tan- 
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gential (parallel and perpendicular to the tube axis). We 
see that the radial modes correspond to frequencies be- 
low 800 cm -1 , corresponding roughly to the out of plane 
bands of graphene. The softening of the frequencies in 
the ab initio calculation are apparent in this figure, spe- 
cially for the higher frequency modes. Also the upper 
limit of radial vibrations is lowered by about 100 cm -1 in 
the ab initio calculation compared with the zone-folding 
results, for the (4,4) tube. 

Several worksEOo have recently made use of the fre- 
quency of the higher optically active modes as an exper- 
imental signature of the tube radii. In graphite there is 
only one Raman active optical mode (with zero wave vec- 
tor) at 1580 cm -1 , which in our calculation appears at 
1690 cm -1 . In the nanotubes, it splits into multiple peaks 
which are originated in the zone folding of the graphene 
bands. The frequencies of these q z =0 tube modes de- 
pend on radius and chirality. In the zone-folding scheme, 
these modes sample the corresponding optical bands of 
graphite, with wave vectors q n — n/r (n — 0,1,2,...) 
along the circumference direction. Of these, only a small 
number (jiadependent of tho-tube radius) are Raman or 
Ir active.BcZl Kasuya et al&B were able to measure the 
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frequency of the highest Raman modes for tubules with 
different radii, and found that these corresponded to the 
n = 1 LO band, and n = 1 and n — 2 TO bands of 
graphite (in order of decreasing frequencies). They found 
the value of the observed frequencies to be in excellent 
agreement with the direct zone-folding results. The ra- 
dius of the narrower tubes in their sample was 5.5 A, 
which, assuming an armchair conformation, would cor- 
respond to the (8,8) tubule. It is, therefore, interesting 
to see whether these results are modified for tubes with 
smaller radii, as was the case for the breathing mode dis- 
cussed above. 

Fig. fj] (a) shows the frequencies of the higher optical 
modes obtained from the ab-initio calculation, whereas 
Fig. fj] (b) shows the results from the zone-folding ap- 
proach. Several facts are worth noticing in the compar- 
ison. As was the case for the breathing mode, the (4,4) 
tube deviates very significantly from the zone-folding be- 
havior. The maximum ab-initio frequency is about 50 
cm -1 lower than the zone- folding prediction, and the 
symmetries and activities of the higher frequency peaks 
are very different. For the larger tubes, it seems that 
the ab-initio results tend to confirm the symmetry as- 
signments of the zone-folding approach, supporting the 
analysis of experimental results in terms of the simple 
zone-folding scheme. The apparent softening of the max- 
imum frequencies for the (10,10) tubule in Fig. [?] (a) is 
due to numerical error (the 10 cm -1 error bar discussed 
above, possibly slightly larger for the largest tube). The 
maximum frequency for this tube should approach the 
one obtained by the zone-folding method since the cur- 
vature is the smallest. This serves as a measure of the 
accuracy of our calculation. 



IV. CONCLUSIONS 

We have presented the results of ab initio calculations 
for single-wall carbon nanotubes with different chiralities 
and radii, addressing structural, elastic, and vibrational 
properties. The cases studied include (n, n) tubes (with 
n ranging from 4 to 10), and the (8,4), and (10,0) tubes. 
We have also presented a detailed comparison with the 
results of other usual approaches like elasticity theory for 
elastic properties or the zone-folding approach for the vi- 
brations. Our results serve to validate most of the pre- 
dictions of these simpler theories, and to point out their 
limits of aplicability. The following conclusions can be 
drawn from our results: 

(i) Relaxation effects due to the tube curvature are 
small in general. The incquivalent bonds in a tube en- 
hance their differences in bond length and angles with de- 
creasing tube radius. The symmetry inequivalent bonds 
give rise to a small shift in the Fermi surface location for 
the (n, n) tubes, mainly related to the lower symmetry of 
the tubes as compared with graphene, and only slightly 
modified by the structural relaxation. 



(ii) The strain energy follows the ^ law expected from 
elasticity theory quite accurately for tubes as narrow as 
(4,4). For armchair tubes, which have slightly lower 
strain energy than other chiralities, the constant has a 
value of a = 2.00 eVA/atom. 

(m) Sensible definitions of the Young modulus are used 
for two different geometries: multi-wall and single- wall 
tubes. In the former case the values are very similar 
to the one of graphite. Single-wall tubes show values 
smaller than graphite. In any case, we propose the elastic 
constant per unit mass as the relevant quantity, since it 
does not depend on the geometry of the system. This is 
shown to be quite similar to the correspondent quantity 
in graphite, for all the studied tubes, and larger than in 
any known fiber. 

(iv) The Poisson ratio also retains graphitic values 
except for a possible slight reduction for small radii. 
It shows a chirality dependence: (n, n) tubes display 
smaller values than (10,0) and (8,4). Our results for the 
(n, n) tubes are consistent with the experimental basal 
Poisson ratio of graphite. 

(v) The phonon bands behave as expected from sim- 
ple schemes, except for slight deviations which become 
more important for the narrower tubes. The zone-folding 
analysis gives a good qualitative (and sometimes quan- 
titative) picture of many of the properties studied here, 
except for known dcfficiencies in the low-frequency vibra- 
tional spectra. For the smallest radii, the zone-folding de- 
scription of the high-frequency vibrations is insufficient, 
too. 

(vi) The breathing mode follows the A/r law predicted 
by graphene-derived force-constants calculations. The 
obtained value of A is consitent with that calculated from 
the in plane stretching elastic constant of graphene. It, 
however, seems to soften with respect to the expectations 
for the smallest radii tubes, like (4,4). A similar softening 
is observed for the twiston modes, whose sound velocity 
diminishes for decreasing radii. 

(yii) The high-frequency optic modes are sensitive to 
the kind of tube and to its radius. The frequencies of the 
highest modes tend to diminish with decreasing radii by 
effect of the curvature. 
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